n_unique=unique(n_soll);
n_min=min(n_unique);
summe=0;
for i=1:length(n_soll)
    if n_soll(i)==n_min
        summe=summe+1;
    end
end
M_auswert=sort(M_mess(1:summe));

[XI YI]=meshgrid(n_unique,M_auswert);
griddata(n_mess,M_mess',eta,XI,YI,'nearest');

M_border_positiv=[0 1.1*max(M_mess)];
M_border_negativ=[0 1.1*min(M_mess)]
for j=1:length(n_unique)
    n_akt=n_unique(j);
    M_max=0;
    M_min=0;
    for k=1:length(n_soll)
        if n_soll(k)==n_akt
            M_akt=M_mess(k);
            if M_akt>M_max
                M_max=M_akt;
            end
            if M_akt<M_min
                M_min=M_akt;
            end
        end
    end
    M_border_positiv=[M_border_positiv;[n_akt M_max]];
    M_border_negativ=[M_border_negativ;[n_akt M_min]];
end
M_border_positiv=[M_border_positiv;[max(n_unique) 1.1*max(M_mess)]];
M_border_negativ=[M_border_negativ;[max(n_unique) 1.1*min(M_mess)]];

figure; hold on; grid off;

linien1=linspace(min(eta),mean(eta)-mean(eta)*0.05,5);
linien2=linspace(mean(eta)+mean(eta)*0.05,max(eta),10);
linien1=linien1*100;
linien2=linien2*100;
linien1=round(linien1);
linien2=round(linien2);
linien1=linien1/100;
linien2=linien2/100;

[C h]=contourf(XI,YI,ans,[linien1 linien2]);
text_handle = clabel(C,h,linien2);
fill(M_border_positiv(:,1),M_border_positiv(:,2),'w') 
fill(M_border_negativ(:,1),M_border_negativ(:,2),'w') 
t=colorbar;
set(get(t,'title'),'String', 'My Title');